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Nature program 


Adrian Owen reveals how to make a Black Spleenwort fern, the 


plant you'll never have to water, and explains what all this has to 
do with the Video Modem 


n January, a C program 

was included in this col- 

umn which will generate 

the Sierpinski triangle. 
This month I’ve included a pro- 
gram that generates a Black Spleen- 
wort fern. Both the triangle and 
the fern are fairly complex images, 
yet the data needed to generate 
them is quite minimal. Look closely 
at the fern program. It picks up 
100,000 random numbers and 
generates 100,000 X,Y coordinates 
from the data stored in arrays A to 
F. The probability array acts as a 
weighting so that certain arrays 
are selected more often than the 
others. The algorithm itself is the 
same as the one used to generate 
the Sierpinski triangle. These ar- 
rays form what is known as an 
‘iterated function set’ or IFS; the 
algorithm which generates these 
images from their data is known 
as the ‘Random Iteration Algo- 
rithm’. 

The Black Spleenwort fern is a 
real plant. It is odd that its rela- 
tively complex shape can be ex- 
pressed in a handful of numbers — 
nature’s numbers, perhaps? Ten 
thousand years ago, all images to 
be found in the world were natural 
images. These images — such as 
trees, mountains and clouds — are 
difficult to generate on a computer 
using conventional computer 
graphic techniques, such as gener- 
ating geometrical shapes, lines and 
spheres. However, any of these 
natural images can be generated 
relatively easily from a set of frac- 
tal numbers. The sets of numbers, 
or iterated function sets, were 
looked at in detail last month. It 
could be said that natural objects 
tend to be made up of copies of 
similar objects to themselves, just 
as any branch on the fractal fern is 
a copy of the whole fern. 

The image of the fern is gradu- 
ally built up by the dancing dot, as 
can be seen in the six fern pictures. 
At first, the image appears as a 
series of random dots, but gradu- 
ally the random points controlled 
by the set of IFS data give the 
deterministic geometrical shape. 
It is important to note that al- 


The Black Spleenwort fern 
is generated using the same 
algorithm as the Sierpinski 
triangle. Notice how the 
image starts off as a 
random collection of dots 
which gradually builds up 
into a structured shape 


though the algorithm works by a 
series of random numbers, the 
resulting image is the same every 
time. The random displacement of 
the dot can be likened to the way 
in which the human eye beholds a 
new object. It darts from point to 
point, building up a picture as it 
goes. 


Man-made images 
On the other hand, man-made 
images which have uniform shapes 


Order from chaos 


and straight lines, such as chim- 
neys, wheels and buildings, can be 
generated using conventional 
computer graphics because the 
squares and spheres that consti- 
tute them can be generated quite 
easily. As the images get more 
complex, so the library of images 
required to generate them becomes 
more complicated — hence the need 
for fractal objects. These complex 
fractal objects, however, can be 
expressed by a handful of num- 


The Sierpinski triangle 
consists of three half-size 
copies of itself, each 
generated by a separate 
code. You can produce 
this image with the pro- 
gram in last month’s 
Advanced Systems 


bers such as the arrays in the Black 
Spleenwort fern program. The 
bitmap for the fern takes up 20,000 
bytes on my disk, but the four 
arrays that make up its IFS take 
100 bytes. 

If Neolithic man had used a 
computer rather than painting in 
the Altamira caves, the natural 
images which were his only sub- 
jects would not have taken up 
much room on his fixed stone 
disks, as they would have 139 
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lent themselves to these elegant 
sets of fractals. 


The Video Modem 

Last month, I began looking at a 
device called the Video Modem. 
Designed by Michael Barnsley and 
Alan Sloan of Iterated Systems in 
Georgia, the device will compress 
video images (256x256 pixels with 
eight bits for colour or greyscale) 
by an average factor of 64:1 and 
transmit 30 of them down a line in 
a second, fast enough to transmit 
motion pictures in real-time down 
a telephone line. It works by break- 
ing the image up into distinct areas 
and then expressing these sub- 
images as a group of iterated 
function sets. The result is that an 
image can be generated from a 
relatively small amount of data, 
compared to the storage needed to 
hold the whole of the bitmapped 
image. 

a We have looked at the task of 
generating an image from its IFS 
in terms of the triangle and fern 
programs. You simply supply the 
iterated function set as a set of 
arrays and the programs, which 
work by the random iteration func- 
tion method, regenerate the im- 
age. The transmitted data, then, is 
the set of iterated function sets 
that constitute the image, and one 
of the functions of the modem is to 
put all sets of incoming iterated 
function sets through a random it- 
eration function in order to regen- 
erate the image and display it at 
the receiving end. 

The other function of the 
modem is to calculate an IFS for 
any given image — ie, to do the 
compression. It enlists the aid of 
the ‘Collage theorem’. 


“he Collage theorem 

cast month, I talked about affine 
transformations. In the fern pro- 
gram, the IFS that generates the 
fern consists of four affine trans- 
formations, the data of the first 
one being [a=0, b=0, c=0, d=0.16, 
e=0, f=0], defining a transforma- 
tion A of any point (X,Y) as 
follows: 


A(x, y) = (ax + by +e, cx + dy +f) 
A(x, y) = (Ox + Oy + 0, Ox + 0.16y +0) 


So during the random iteration 
function, say a point (10, 20) arises, 
and the first affine transformation 
is chosen: 


A(x, y) = (0*10 + 0*20 + 0, 0*10 4 
0.16*20 +0) 


which gives the next point as (0, 
3.2). 


Watchmen 


Here’s an interesting affine trans- 


formation: 


B(x, y) = (1/2x + 1/4y +1, 1/4x + 1/2y 
+2) 


If this transformation is applied to 
every point on the Watchmen 
badge, the resulting image appears 
flattened and on its side, ie con- 
tracted. This is the meat and bones 
of the Collage theorem. By choos- 
ing a good approximation of the 
shape to generate, it is then pos- 
sible to calculate the affine trans- 
formation that will contract the 
approximated shape into the re- 
quired one. 

In the Watchmen badge ex- 
ample, mark three points on the 


standard badge and three points 
on the contracted badge, then 
calculate the affine transforma- 
tion which will generate this new 
image from the old by solving 
these linear equations: 


a, bande 
ala+a2b+e=a1 (1) 
bla+b2b+e=b1 (2) 
cla+c2b+e=c (3) 
c, dandf 


aic+a2d+f=a2 (1) 


BLACK SPLEENWORT FERN 


#include <math.h> 
#include <graphics.h> 
#include <stdio.h> 
#include <conio.h> 
#include <dos.h> 


int graphdriver=DETECT, graphmode, lw, op,xc,yc,ang; 
intj,k,xscale, yscale,xoffset, yotfset, pr,p[4],pk[4]; 


long unsigned int i; 


float a[4], b[4], c[4], d[4], e[4], f[4], x, y, newx; 


main() 
{ 


int px, py; 


/* Initialise the graphics */ 


initgraph( &graphdriver, &graphmode,””); 


/* Plug the data into the Iterated Function Set */ 
a[0] = 0; a[1] = .20; a[2] = -.15; a[3] = .85; 

b[0] = 0; b[1] = -.26; b[2] = .28; b[3] = .04; 

c[0] = 0; c[1] = .23; c[2] = .26; [3] = -.04; 

d[0] = .16; d[1] = .22; d[2] = .24; d[3] =.85; 

e[0] = 0; e[1] = 0; e[2] = 0; e[3] = 0; 

{[0] = 0; f[1] = 1.60; f[2] = .44; f[3] = 1.6; 

pl0] = 328; p[1] = 2621; p[2] = 4915; p[3] = 32767; 


/* Set the controls for the heart of the sun */ 


xscale = 35; 
yscale = 35; 
xoffset = 300; 
yoffset =0; 


x=0; 
y=0; 


for(i=1; i<=100000; i++) 
{ 


j = rand(); 


k = (j<p[0]) ? 0: ((j<p[1]) ? 1 : ((i<p{2]) ? 2 : 3)); 
newx = (a[k] * x+b[k] * y+e[k]); 
y = (c[k] * x+d[k] * y+f{k]); 


xX = newx; 


px = x* xscale + xoffset; 
py = (350 - y * yscale + yoffset); 


if((px>=0) && (px<640) && (py>=0) && (py<350)) 
putpixel(px, py, 7); 


b1c+b2d+f=b2 (2) 
cle+c2d+f=c2 (3) 


Voila! Given any shape we wish to 
generate, find a good approxima- 
tion to it in the library of shapes, 
and calculate the affine transfor- 
mation to change the library shape 
into the required one. 

But what if there isn’t a good 
approximation? In this case, the 
task is a little more difficult, but 
here’s where the ‘collage’ comes 
into play. A collage of known 
shapes is assembled so that the 
union of all of them gives a good 
approximation of the image we 
wish to generate. Calculate the 
affine transformations that will 
produce each of the new shapes 
from each of the old and what 
remains is a number of affine 
transformations, the union of 
which results in the new image 
being generated. This group of 
affine transformations is the iter- 
ated function set for that image, 
which will generate it when put 
through the random iteration func- 
tion. 


The library 

The missing piece of this jigsaw is 
how the Video Modem spots the 
best collage to use — in other 
words, how it finds the best ap- 
proximation of shapes from the 
library to make up the required 
image. 

In January, I explained that 
the modem first preprocesses a 
video image to break it up into a 
number of sub images. The task is 
then to find the iterated function 
sets for each of these sub images. 
But the sub images themselves 
could potentially be made up of a 
collage of any number of images. 
How can the modem tell if four or 
54 pieces are required to give the 
best approximation, and how 
exactly does the matching of con- 
tractive library images to the re- 
quired bit of a collage take place? 
Conventionally, it is possible to 
envisage a point-and-drag arrange- 
ment using a mouse to drag a 
suggested library image until it fits 
over the required image. This is, of 
course, manual, so how does the 
Video Modem manage to do it in 
real-time? Mr Barnsley and Mr 
Sloan are obviously keeping that 
one under their hats. 


Adrian Owen is currently 
working on an iterated 
function set to generate a 
fully-stocked bar 
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Olivetti PC Amstrad PC 
MONO COLOUR 
PCS 86 SD 549 735 
PCS 86 DD 649 825 ioeoe 
PCS 86 HD 825 1009 


MONO COLOUR 
349 429 
1512 DD 429 505 


1640 SD 440 530 
PCS 286 DD 1099 1640 DD 499 582 
PCS 286 SD2 1 196 1640 HD30 675 759 836 
PCS 286 SD40HD 12” VGACOL 14” VGACOL 
R 687 777 


All olivetti PCs come g 9 2086 SD 519 599 
site maintenanc 2086 DD 585 669 749 837 


Frye! 86 HD3O 745 839 923 994 
PRINTERS 


PHILIPS NMS1432 

STAR LC-10 MARK II 
STAR LC-10 

STAR LC-10 COLOUR 
STAR LC24-10 

STAR XB24-10 COLOUR 
PANASONIC KXP11-24 


HP DESKJET ull-ra d Philips PCs come 

HP PAINTJET with 12 months on- 
CITIZEN 120D Ms site maintenance 

EPSON LX-400 CORD 

EPSON LX-850 D 

EPSON LQ-500 249 9115 20HD 


OKIMATE 20 159 9115 20HD ECD 
LASERS AT SERIES 
PHILIPS LASER 1299 9125 20HD SD MD 
CANON LASER MARK II 1495 9125 20HD SD CD 
PANASONIC KX-P4450 1395 9125 HD SD ECD 
HP LASERJET II 1459 
NEC POSTSCRIPT 2589 ] 
NEC POSTSCRIPT ses MR DIAMONDs SPECIAL! 


OLIVETTI PG-303 


Twin drive CGA - £499 + VAT 


32 Mb Hard Card - £199 
Exclusive mobile rack — the ultimate 
security! Carry your hard disk with you! 
£69.95 inc vat 


PC3 DD MONO 

PC3 DD EGA 

PC3 DD 32HD MONO 
PC3 DD 32HD EGA 899 
PC4 BASE ONLY 1199 
PC4 MONO 1299 


PC4 EGA COLOUR 1449 
PC4 VGA COLOUR 1499 


ae onsite maintenance. - A BRING THIS AD TO A 
951 
Atari Portfolio — £189 D’GPT, ONE FOR £4.95 


Accessories: 
Serial adaptor £59.95 
Parallel adaptor £39.95 
Memory upgrade to 640K — £179.95 
32K memory card — £49.95 
64K memory card — £89.95 
128K memory card - £129.95 
Card drive including PC read/write —- £69.95 
Accessory prices include vat 


